Entanglement Spectrum as a Generalization of Entanglement Entropy: Identification 
of Topological Order in Non-Abelian Fractional Quantum Hall Effect States 
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We study the "entanglement spectrum" (a presentation of the Schmidt decomposition analogous 
to a set of "energy levels" ) of a many-body state, and compare the Moore- Read model wavefunction 
for the v = 5/2 fractional quantum Hall state with a generic 5/2 state obtained by finite-size 
diagonalization of the second-Landau-level-projected Coulomb interactions. Their spectra share a 
common "gapless" structure, related to conformal field theory. In the model state, these are the 
only levels, while in the "generic" case, they are separated from the rest of the spectrum by a clear 
"entanglement gap", which appears to remain finite in the thermodynamic limit. We propose that 
the low-lying entanglement spectrum can be used as a "fingerprint" to identify topological order. 
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There has been increasing recent interest in using 
quantum entanglement as a probe to detect topological 
properties of many-body quantum states [l], [2, HI, |4j, in 
particular the states exhibiting the fractional quantum 
Hall effect (FQHE) @, 0, 0, H, M, Efll, El El, El ■ Among 
various measures of quantum entanglement, the entan- 
glement entropy has by far been the favorite H, lil, 3 • 
By partitioning a many-body quantum system into two 
blocks, the entanglement entropy is defined as the von 
Neumann entropy of the reduced density-matrix of either 
one of the two blocks, and is a single number that can 
be obtained from knowledge of the density-matrix eigen- 
value spectrum. The density matrix may be written in 
the form p = exp(— H), so that the entanglement entropy 
is equivalent to the thermodynamic entropy of a system 
described by a hermitian "Hamiltonian" H at "temper- 
ature" T = 1; in the case of a weak entanglement, the 
"excited states" eigenvalues of H are separated from the 
ground state eigenvalue by a large "energy gap" that be- 
comes infinite in the limit of a simple product state with 
vanishing entanglement entropy. In this Letter, we point 
out that the spectrum of the "Hamiltonian" H, which 
we call the "entanglement spectrum", reveals much more 
complete information than the entanglement entropy, a 
single number. 

We examined the entanglement spectrum of the ideal 
Moore- Read (MR) state with Landau-level filling fraction 
v = 1/2, and found that its structure appears to coincide 
with that of the spectrum of the associated nonabelian 
conformal field theory (CFT). In particular, the spec- 
trum is "gapless" , and (up to a limit determined by the 
compactification of the system to have a finite area) has 
the same count of states (characters) at each momentum 
(Virasoro level), although there is some splitting of the 
"energies" . The MR entanglement spectrum is also "non- 
generic" , as it contains far fewer levels than the number 
expected from counting Hilbert-space dimensions; this is 
a common property of model wavefunctions, which lack 
part of the "zero-point fluctuation" of a generic state. 

We also examined the "realistic" i/ = 5/2 = 2+l/2 
state obtained by diagonalizing the Coulomb interaction 



projected into the half-filled second Landau level, and 
observed that all the extra levels expected in a generic 
entanglement spectrum are present, but are separated by 
a finite gap from "gapless" CFT-like low-lying states with 
essentially the same structure as those of the Moore- Read 
state. This gap appears to remain finite in the thermody- 
namic limit. Thus the low- lying entanglement spectrum 
is essentially a "fingerprint" which allows the associated 
CFT (which characterizes the topological order) to be 
identified. 

If Landau-level mixing is ignored, the v — 5/2 FQHE 
system is equivalent to a v = 1/2 lowest Landau level 
with interaction pseudo-potentials [lH corresponding to 
a simple Coulomb interaction projected into the second 
Landau level (we did not refine this to include the quan- 
tum well form-factor, which may be important for quanti- 
tative modeling) . In spherical geometry ll| , the number 
of electrons (N e ) and the total number of Landau level 
orbitals (N or t,) are related by N or b — 2N e — 2. In second- 
quantization, a many-body state can be written in the 
basis of orbital occupations; we used the L z eigenstate 
basis to divide the spherical surface at a line of latitude 
into two regions, so the N or i, orbitals are partitioned into 
^orb around the north pole, and N^ rb around the south 
pole (this is spatially the sharpest cut consistent with 
projection into a Landau level), which is a partition of 
the Fock space TL into two parts Ha ®~Hb- 

A Schmidt decomposition (equivalent to the singular 
value decomposition of a matrix) of a many-body state 



gives 
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where exp(-|&) > 0, \vj a ) € Ha, \i> B ) e Kb, and 

{^aWa) = {^bWb). = S ivm g exp(-±&) as the sin- 
gular values and |^) and \ip B ) the singular vectors. If 
the state is normalized, X)i ex P(~ &) = 1, but ^ is n °t 
necessary to impose the normalization condition. 

The £,'s are "energy levels" of a system with ther- 
modynamic entropy at "temperature" T = 1 equivalent 
to the entanglement entropy, S — ex P ( — CO, which 
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FIG. 1: The complete entanglement spectra of the JV e = 16 and N or t = 30 Moore- Read state (only the relative values of £ and 
Lf are meaningful). 



has been shown to contain information on the topolog- 
ical properties of the many-body state [12j. The full 
structure of the "entanglement spectrum" (logarithmic 
Schmidt spectrum) of levels £j contains much more infor- 
mation than the entanglement entropy S, a single num- 
ber. This is analogous to the extra information about 
a condensed matter system given by its low-energy ex- 
citation spectrum rather than just by its ground state 
energy. 

Because the FQHE ground state is translationally and 
rotationally invariant (with quantum number L to t = on 
the sphere) , and the partitioning of Landau-level orbitals 
conserves both gauge symmetry and rotational symme- 
try along the z-direction, in either block A or B both the 
electron number (Nf and Nj?) and the total z-angular 
momentum (Lf and Lf ) are good quantum numbers 
constrained by Nf + N? = N e , Lf + Lf = 0. The en- 
tanglement spectrum splits into distinct sectors labeled 
by Nf and Lf. 

In the thermodynamic limit, the MR model state can 
be represented by its "root configuration" 15j, which 
has occupation numbers "11001100 ••• 110011" , with a 
repeated sequence •••1100---; in spherical geometry, 
this is terminated by "11" at both ends. This is the 
highest-density "MR root configuration", which we de- 
fine as an occupation-number configuration satisfying a 
"generalized Pauli principle" that no group of 4 consec- 
utive orbitals contains more than 2 particles (this rule 
also applies to MR states with quasiholes, and gener- 



TABLE I: The numbers in the parenthesis are values of 
(N^. b , Nf), respectively for each system and partitioning as 
specified. 



N e 


P[0|0] 


P[0|1] 


P[l|l] 


10 


(7,4) 


(8,4) 


(9,5) 


12 or 14 


(11,6) 


(12,6) 


(13,7) 


16 


(15,8) 


(16,8) 


(17,9) 



ates the CFT edge spectrum of a finite MR droplet on 
the open plane 15|.) When the MR state is expanded 
in the occupation-number basis, the only configurations 
(Slater determinants) present are those obtained by start- 
ing from the root configuration, and "squeezing" pairs 
of particles with L z — mi,rri2 closer together, reducing 
| mi — jtt-2 | , while preserving mi + mi [la ]. 

From the root configuration, we see that there are three 
distinct ways of partitioning the orbitals: (i) between two 
0's; (ii) between two l's; or (iii) between and 1 (the par- 
titioning between 1 and is equivalent to that between 
and 1 by reflection symmetry). We use symbols P[0|0], 
P[l|l], and P[0|1] to represent the three cases, respec- 
tively. This will correspond to choosing one of the three 
sectors of the associated conformal field theory. For finite 
systems, we always try to draw the boundary of the par- 
titioning cither on the equator (if possible) , or closest to 
the equator but in the southern hemisphere. Moreover, 
we can associate a "natural" value to for a particu- 
lar partitioning, i.e., the total number of l's in the root 
occupation sequence on the left-hand-side to the bound- 
ary. In this Letter, it is sufficient to consider only levels 
whose Nf is exactly this natural value. Table [I] describes 
the precise meaning of these symbols for systems that are 
considered here. 

Figure (p} shows the spectra for each of the three dif- 
ferent ways of partitioning, for the Moore-Read state at 
N e = 16 and N or b = 30. The spectrum not only has 
far fewer levels than expected for a generic wavefunction, 
but also exhibits an intriguing level-counting structure 
(as a function of Lf and Nf-) that resembles that of the 
associated conformal field theory of the edge excitations. 
Intuitively, this is because the boundary of the Landau- 
level partitioning indeed defines an edge shared by region 
A and B. 

In the intuitive picture, the quantum entanglement 
between A and B arises from correlated quasihole exci- 
tations across the boundary along which the partition- 
ing is carried out. Any quasihole excitation in region 
A necessarily pushes electrons into region B, and vice 
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FIG. 2: The low-lying entanglement spectra of the N e = 16 and N or t, — 30 ground state of the Coulomb interaction projected 
into the second Landau level (there are levels beyond the regions shown here, but they are not of interest to us). The insets 
show the low-lying parts of the spectra of the Moore- Read state, for comparison [see Figure {TJ]. Note that the structure of 
the low-lying spectrum is essentially identical to that of the ideal Moore- Read state. 



versa. However, the electron density anywhere on the 
sphere must remain constant, which can be achieved if 
the quasihole excitations in A and B are correlated (en- 
tangled) . This gives the empirical rules of counting the 



levels. Take the spectrum in Figure ( 1(a) ) as an example 



The partitioning P[0|0] results in the root configuration 
110011001100110 on the northern hemisphere (region A), 
and it corresponds to the single "level" at the highest pos- 
sible value of Lf = Lf max = 64. We measure the Lf by 



its deviation from Lf max , i.e. AL := Lf >max -Lf, which 
has the physical meaning of being the total z-angular 
momentum carried by the quasiholes. At AL — 1, the 
levels correspond to edge excitations upon the AL = 
root configuration. There is exactly one edge mode 
in this case, represented by the MR root configuration 
110011001100101. 

The number of AL = 2 levels can be counted in exactly 
the same way. There are three of them, of which the root 
configurations are 

1100110011001001 
1100110011000110 
1100110010101010 

while for AL — 3, the five root configurations are 

11001100110010001 
11001100110001010 
11001100101001100 
11001100101010010 
11001010101010100 

The counting for the levels at small AL for P[0|1] and 
P[l|l] can be obtained similarly. 

For an infinite system in the thermodynamic limit, the 
above idea gives an empirical counting rule of the number 
of levels at any AL, i.e., it is the number of independent 
quasihole excitations upon the semi-infinite root config- 
uration uniquely defined by the partitioning. For a finite 



system, this rule explains the counting only for small 
AL; for large AL, the finite size limits the maximal an- 
gular momentum that can be carried by an individual 
quasihole. Therefore the number of levels at large AL 
in a finite system will be smaller than the number ex- 
pected in an infinite system. Not only is this empirical 
rule consistent with all our numerical calculation, but 
it also explains why P[0|0] and -P[0|1] have essentially 
identical low-lying structures. This is because the (semi- 
infinite) configuration "• • • 1100110" is essentially equiva- 
lent to "• • ■ 11001100" (with an extra "0" attached to the 
right). We expect that P[0|0] and P[0|1] become exactly 
identical in the thermodynamic limit. 

For completeness, we list the root configurations asso- 



ciated with the first few low- lying levels in Figure (1(c)) 



AL = : 11001100110011001 
AL = 1 : 110011001100110001 
110011001100101010 
AL = 2 : 1100110011001100001 
1100110011001010010 
1100110011001001100 
1100110010101010100 

Figure ^ shows the spectra of the system of the same 
size as in Figure (JTJ), i.e., N e = 16 and N or b = 30, but for 
the ground state of the Coulomb interaction projected 
into the second Landau level, obtained by direct diag- 
onalization. Interestingly, the low-lying levels have the 
same counting structure as the corresponding Moore- 
Read case. We identify these low-lying levels as the 
"CFT" part of the spectrum, in contrast to the other 
generic, non-CFT levels that are expected for generic 
many-body states. At relatively small AL (up to a limit 
which grows with the size of the system), the CFT lev- 
els are separated from the generic levels by a clear gap, 
which we define as the distance from the average of the 
CFT levels to the bottom of the generic levels. 
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FIG. 3: Entanglement gap as a function of l/N. So is the gap at AL — 0, i.e., the distance from the single CFT level at 
AL = to the bottom of the generic (non-CFT) levels at AL = 0. At AL — 1, 2, the gap <5i,2 is defined as the distance from 
the average of the CFT levels to the bottom of the generic levels. See Table U for the details of various partitionings. 



Figure §5§ shows the value of this gap (at AL = 0,1,2 
respectively) as a function of the size of the system, based 
on which we speculate that the gap between the CFT 
and non-CFT levels remains finite in the thermodynamic 
limit for all AL. The observed fact that the structure of 
the low-lying spectrum is essentially identical to that of 
the Moore- Read state, as well as the existence of the "en- 
tanglement gap", serve as evidence that the "realistic" 
v = 5/2 FQH states is indeed modeled by the Moore- 
Read state. 

Assuming that the gap does remain finite in the ther- 
modynamic limit, characterization of the entanglement 
spectrum is a reliable way to identify a topologically 
ordered state. While finite-size numerical studies of- 
ten show impressive (e.g., 99%) overlaps between model 
wave- functions (Laughlin, Moore-Read, etc.) and "re- 
alistic" states at intermediate system sizes, this cannot 
persist in the thermodynamic limit. Furthermore, the 
entanglement spectrum is a property of the ground state 
wave- function itself, as oppose to the physical excitations 
of a system with boundaries, so allows direct comparison 
between model states and physical ones. 

The asymptotic behavior of the characters of a CFT 
(the count of independent states at each Virasoro level, in 
each sector) defines both the effective conformal anomaly 
c of the CFT, and the quantum dimension of each sec- 
tor. To the extent that there is a clear separation of the 
"gapless" , low-lying CFT-like modes and the generic (but 
"gapped" ) modes, one can count the number of "gapless" 
modes as a function of momentum parallel to the bound- 
ary separating the two regions. For a finite-size system, 
these will match the CFT characters up to some limit 
that grows with system size. These numbers are integers, 
so are not subject to numerical error, and in principle, 
both c and the quantum dimensions can be extracted 
from their behavior as the system size grows. 

As a critical point is approached, the "entanglement 



gap" may still be finite but well below the "temperature" 
T = 1 at which the von Neumann entropy is evaluated. 
We suggest that the direct study of the low-lying entan- 
glement spectrum is a far more meaningful way to char- 
acterize bipartite entanglement. Equivalently, the T — > 
"low-temperature" entropy of the modified family of den- 
sity matrices p^ 1 ^ T ^ may prove useful, as this corresponds 
to the thermodynamic entropy of the entanglement spec- 
trum at temperature T. 

We thank B. A. Bernevig for valuable comments and 
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MRSEC DMR02-13706. 
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